Decoherence-free quantum information in Markovian systems 
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Decoherence in Markovian systems can result indirectly from the action of a system Hamiltonian 
which is usually fixed and unavoidable. Here, we show that in general in Markovian systems, 
because of the system Hamiltonian, quantum information decoheres. We give conditions for the 
system Hamiltonian that must be satisfied if coherence is to be preserved. Finally, we show how to 
construct robust subspaces for quantum information processing. 



Decoherence remains the most important obstacle to 
experimental realizations of quantum processors. One 
well-developed method of counteracting the effects of de- 
coherence is to encode quantum information (QI) into 
decoherence-free subspaces and subsystems (DFSs) [l|, 0, 
& 0, H, @, • This form of passive error correction has 
been well-studied, and has led to (approximate) realiza- 
tions of DFSs. For example, in ion traps Kielpinski et 
al. used a decoherence-free (DF) state of two trapped ions 
to enable encoded information to be stored longer than 
its unencoded counterpart Q, and in an optical system 
Kwiat et al. prepared a similar DF state using paramet- 
ric down-conversion 0. These proof-of-principle experi- 
ments have shown that encoding QI into DFSs improves 
storage lifetimes, and have partly justified the extensive 
theoretical investigations into DFSs. 

The (strict) requirement for infinite-lifetime (DF) 
quantum information storage is that all qubits must be 
symmetrically coupled to the environment [3, 0, 0, 0, 0, 0, 
[7[. Most theoretical results regarding DFSs rely on this 
assumption, which is possible to obtain only for qubits 
that are colocated. As the number of qubits increases 
this approximation becomes less tenable. Furthermore, 
it has recently been shown that infinite-lifetime quantum 
information storage is not possible if the physical qubits 
of interest are not colocated [13]. So, at best we are 
left with the regime in which all physical qubits approxi- 
mately experience the same Markovian environmcnt-the 
qubits are very close together, but not colocated. A full 
characterization of decoherence in this regime is the pur- 
pose of this paper. 

So, we examine the hitherto neglected case of decoher- 
ence in physical systems for which the exact symmet- 
ric coupling between the qubits is pcrturbatively bro- 
ken. One physical example of this is given by closely 
spaced dipole-coupled qubits that approximately satisfy 
the requirements for Dicke superradiance [Til ] . For a full 
analysis of this experimentally accessible regime within 
the Lindblad master equation, the effects of both the 
Lamb-shift-type Hamiltonian and the decoherence op- 
erators must be included. The unitary evolution gen- 
erated by the Lamb-shift-type Hamiltonian can cause 
DFQI to evolve into non-DF states, and so decay via 



the action of the decoherence operators. Another exam- 
ple is a Hciscnbcrg-typc interaction in a spin lattice. We 
make no assumptions with regards to the Lamb-shift- 
type Hamiltonian, and find that for Markovian systems 
stable quantum information is rare. We relax the re- 
quirement for infinite-lifetime information storage, and 
derive expressions for fidelities which depend on the rel- 
ative strengths of the evolution operators. Finally, we 
show how to construct robust subspaces for quantum in- 
formation processing. 

For the system density-matrix p in the Hilbert space 
7i, the most general description of Markovian dynamics 
for initial decoupling between the system and the bath is 
given by the Lindblad master equation 

p= -i[H,p} + L D (p) (la) 

L d(p) = \ E MIS*, PS/] + [S iP , S/]), (lb) 

where (a^) is a time- independent Hermitian coefficient 
matrix, and Ld is a completely positive map. The 
presence of the decoherence operators Si in Eq. (jlbj) 
means that the evolution may not be unitary. But, if 
L D (p(t)) = then p(t) = e~ iHt p(0)e lHt , so one defines 
DF dynamics as satisfying L_n(p) = 0. One can restrict 
to pure states p = \ip) [5[ to give a sufficient con- 
dition for L D (p) = as S l \ip) = a \ip). Thus, a DFS 
is defined as an invariant subspace M c TL such that 
Si = ct V t/j e M. This condition is guaran- 
teed if M is an eigenspace of Sj with eigenvalue zero, 
which in many cases of interest is true for all genera- 
tors. We note that in two important cases this is the 
only possibility: i) if the Lie algebra generated by Si is 
scmisimple and ii) if all the generators satisfy S™ = 
for some n (are nilpotent). More precisely, a subspace 
A e TL is DF if Si(A) = V i. Although, this condi- 
tion seems somewhat stronger than the usual condition 
of a common invariant eigenspace with eigenvalue Ci (not 
necessarily zero) of Si the analysis below can be gen- 
eralized by transforming the operators Si — > Si — c,7. 
The evolution equation ([1]) has a unitary and nonunitary 
part, so even if p(to) satisfies Lrj(p(t )) = it does not 
ensure that L^(p(t)) = 0, for t > to. This motivates 
the following definition. A DFS A of TL is to be consid- 
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ered complctcly-decohcrence-free (CDF) if for any state 
| a) G A, p(0) = | a) (a | gives L D {p{t)) = 0. This condi- 
tion is weaker than that derived in Ref . [l2[ , but stronger 
than that derived in Ref. 

We illustrate the regime of interest to this paper using 
the example of dipole-coupled qubits. For these qubits, 
there exists a regime for which H = V- . A^S^S 3 . 
This occurs when the qubits are closely spaced [111 ]. 
It happens that for closely-spaced dipole-coupled qubits 
the spatially-dependent interatomic spontaneous emis- 
sion described by the matrix (a^ ) no longer depends on 
index i,j, but the spatially-dependent interatomic coher- 
ent dipole-dipolc interaction docs. This means that there 
exists exact DF states that are acted upon by H. caus- 
ing transitions to non-DF states. This is the case for any 
separation greater than zero. 

We begin with a criteria for CDF dynamics. 

Proposition 1 Let V be a DFS in H: S; -V = V i. A 
necessary and sufficient condition that V contain a CDFS 
W is H • W C W. In particular, H can be diagonalized 
in W. 

Proof Define p'(t) = e iHt p(t)e~ im . The equa- 
tion satisfied by p' is p' = L' D [p'}, where L^fS'J = 
L D [e mt S l e- mt ]. Hence, in this picture p'(0) is DF iff 
L' D [p'(0)] = 0. The generic DFSs are spanned by vectors 
\x) such that \x) = e 4H *S 4 e- lHt \x) = 0. Let W be 
the subspace consisting of all such vectors. This must be 
satisfied for all t, so we have H • W C W Vi. Conversely, 
if this is satisfied then Sie~ lIit \x) = Vi. 

A useful consequence of the proposition is the follow- 
ing. 

Corollary A subspace W is CDF iff (ad(H)) n (Si) • 
W = 0, Vn and i. This is equivalent to the condition 
[H n ,Si] -W = 0. 

Proof The first condition follows from the proposition 
and the identity e 4tadH (S l ) = e 4tH S 4 e- ltH . The second 
condition is proved by induction [l3j |. 

The Pauli matrices {<r x , a y , a z } generate the Lie al- 
gebra su(2) and form a basis (with identity matrix /) 
for the space of observables for each qubit. Let S l a = 
I ® • • • ® <J a ® I ■ ■ ■ ® I, a £ {x,y, z}, with <r a only at 
the i th place. Let S l ± = ± iS», S± = £ 4 S 4 ± and 
S 2 = 2 S l z ■ These operators define a representation 
of su(2) on the system Hilbert space H = ® N C 2 = C 2N , 
and satisfy [S±,SJ = -2S Z and [S+, S_] = S*. We 
denote a tensor basis for C 2N as (8> • ■ ■ |*jv) where 
ik S {0,1} = [ix • • - *iv) 3 an d we write \ji ■ ■ ■ jk) for a 
vector with ones at j\ ■ ■ ■ jk and zeros elsewhere. The 
Hamiltonian is written H = ^ijS l + S J _, where (A^) 
is a Hermitian matrix. Without loss of generality, we 
take A, for A = (Ay), to be real symmetric. We con- 
sider the case where a^- = a, so there is only one Lind- 
blad generator S_. In this instance, Eq. (jlbj) takes the 
form L D (p) = k(S_ /9 S+ - i(S+S_p + pS+S_)). Each 
irreducible representation (irrcp) of TL is generated by a 
unique lowest weight vector with weight — (r— 1) for irrep 



dimension r that satisfies S_ \a) = 0. Note that H leaves 
the weight spaces invariant, so we only consider DFSs of 
fixed weight. 

The subspace V\ is generated by the basis B = {\i) = 
|0 ■ ■ • 010 • • ■ 0} | for 1 in the i th place}. Here, we let H 
also stand for its restriction to V\. Then, H \i) = 
J2k^ki \k). A state |x) = J2i x i I*) i s a lowest weight 
state iff Y^ii a i = 0. There are N — 1 such independent 
vectors which generate a DFS, called here D\. We wish 
to find out whether there are any CDFS ^ inside D\. 
Note that we refer to the CDFS C of D\, which denotes 
the maximal subspace that is the sum of all CDFSs. The 
condition for CDF dynamics is H • C C C . So, for a 
nonzero subspace C to exist it is necessary and sufficient 
that H have an eigenvector in V\. In the fixed basis B, 
we represent an arbitrary |x) = X^i 2 -* I*) as a column 
vector x T = (x±, ■ ■ ■ ,Xn) in C^. It is clear that |x) is 
an eigenvector of H iff x is an eigenvector of A with the 
same eigenvalue. Hence, there will be a nonzero CDFS 
iff A has an eigenvector x such that Tr(x) = where 
Tr(x) = J2i x i- Suppose A has a degenerate eigen- 
value c. Then there are at least two independent vectors 
y and z. If we have Tr(y) = Tr(z) = 0, then there is a 
nonzero CDFS containing at least |y) and |z). Other- 
wise, suppose Tr(z) = k ^ 0. Then, x = y Tr j^ z has 
trace zero and we have a nontrivial CDFS. This gives: 

Proposition 2 A sufficient condition for the existence 
of nonzero CDFS is that the matrix A has a degenerate 
eigenvalue c. If c is m-fold degenerate, then the dimen- 
sion of the CDFS > m - 1. 

Although the main result of this paper concerns the 
rarity of CDFSs, there are two examples of physical sys- 
tems that satisfy CDF dynamics. First, consider four 
qubits in a spatially symmetric configuration, e.g., the 
corners of a square lattice. Then, from symmetry con- 
siderations it is clear that the row (or column) sums of 
the matrix A is constant. This automatically guaran- 
tees the condition in the proposition. There is at least 
one CDFS containing Ylt=i Second, consider two 
dipole-coupled two-level atoms. In the Dicke limit, the 
single-excitation antisymmetric state is CDF because for 
the special case of two atoms, the dipole-dipole inter- 
action does not cause information to leak from the DF 
state to the non-DF state. Consider now the general 
case for a single-excitation. Assume all the eigenvalues 
of A are non-degenerate. We seek a condition on A 
that will ensure the existence of an eigenvector in D\. 
Suppose x = (xi, ■ ■ ■ ,xn) t is such an eigenvector with 
eigenvalue c. Then the corresponding eigenspace is 1- 
dimensional, and we have the following set of equations 
AnXi+- • -+AiAr.Tjv = cxi; AjviXi+> • -+A nn xn = cx N 
; X\ + • • • + xjs; = 0. If we set one of the components, 
say xn = 1, these equations have a unique solution. 
We rewrite the first N equations as (An — c]x\ + ■ ■ ■ + 
Ai.Ar_i.Tjv_i = — Aitv; Ajv_x,ia;i + • ■ ■ + (Ajv-i,jv-i — 
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c)x N -i = -A N -i,n; A N iXi + ■ 
r(c) denote the (N - 1) X (N 



+ Ann — c = 0. Let 
1) matrix such that 



action of H on V m is more complicated: 



cSij, 1 < i, j < N — 1 (for (5^ the Kronecker 



r = A 

(5). The uniqueness of x' = (xi, . . . , XAr_i) T implies that 
r is invertible. Writing d = — (Auv, ■ ■ ■ , An-i,n) T we 
have x' = r _1 (c)d. The last equation Ajvia;i + ■ • • + 
Ann — c = can be written as d T x' — c + Ann- 
Hence, we have d T r _1 (c)d — c + Ann = 0. That is, 
/(c) = d T adj(r(c))d - det (r(c))(c - A NN ) = , where 
adj(r) is such that adj(r)T = det (T)I. From the condi- 
tion Tr(x) = and i« = 1 wc obtain Tr(r _1 (c)d) = -1. 
Hence, g(c) = Tr(adj(r(c))d) + det (T(c)) = 0. So, we 
get two polynomial equations in c whose coefficients are 
functions of An . For a solution to exist the resultant 



of the two polynomials must vanish 14|, and we get a 
polynomial relation i£(Ay) among the A^ which does 
not vanish identically. There are N(N + l)/2 = M in- 
dependent parameters characterizing any real symmetric 
matrix. Hence, the space of all such matrices may be 
identified with M. M . We have just seen that for A to have 
a nontrivial CDFS, it must satisfy (at least) one polyno- 
mial equation. So, there are Hamiltonians with matrix 
A which do not have any CDFS in V\ . We will see that 
this is the norm rather that the exception — Hamiltonians 
with CDFS are rare. 



We illustrate the above analysis with an example. The 
first excited subspace V\ is spanned by{|l) = |001),|2) = 
|010) , |3) = 1 100) } . Let the matrix A corresponding to 
H restricted to V\ be given by A i2 = A 2 i = x$, A13 = 
A31 = X2 and A23 = A32 = x\ with diagonal elements 
zero. The CDFS condition translates to X3 = x\ or X2 
corresponding to eigenvalues — X2 and — x\ respectively. 
The third possibility x\ = X2 gives eigenvalue — X3. 
Fixing X3 = X\ = a the eigenvalues of A are 1 and 
X2 ± \J x\ + 8a 2 . The only possibility for a degenerate 
eigenvalue is when x± — X2 = X3. We conclude that H 
restricted to the first excited subspace V\ will have an 
eigenvector in D\ (DF subspace) iff at least two of its 
off-diagonal entries arc equal, and it will have two eigen- 
vectors in D\ if all three are equal. The parameter space 
of H can be identified with R 3 as x\, X2 and X3 the range 
over the real numbers. Then, the only Hamiltonians with 
eigenvectors in D\ arc characterized by the parameters 
that lie in the planes X\ = X2, x% = X3OT x 2 = X3. So, we 
conclude that the Hamiltonians which leave some DFS 
state in V\ invariant is a negligible fraction of all the 
possible Hamiltonians — in general, quantum information 
will decay. 



We know that the DFS in the m th excited subspace 
V m is spanned by vectors |x m ) = Y^ x ii,-,i m Ki; • ■ ■ 
such that M' = ( m _ 1 ) equations 



N 



(*) — with all indices except i r fixed — are satisfied. The 



H|x m )= ^ x h,-,j m A kn \ki2i3 ■ ■ ■ i m ) + 

ii,...,i m k 

■■■ + J2 A kim 1*1*3*3 •••*!»■ (2) 



We require that |x m ) satisfy H |x m ) = A |x m ). We write 
the matrix representing H restricted to the subspace V m 
as A( m ). The action of H on V m is equivalent to that of 
A( m ) on 1R M for M = ( ) whose coordinates are given 
by ^ii,...,v,...,i m - We write the eigenvalue equations as in 
the previous section. Since the eigenvector must satisfy 
M' equations (*) we use them to write the last M 1 com- 
ponents of such a vector in terms of the first M — M' , 
and substitute in the eigenvalue equation of A^™* 1 . The 
resulting system of equations in M — M' variables must 
have rank less than M — M' for a nontrivial solution 
to exist. Let A'( m ) be the square matrix of the coef- 
ficients of the first M — M' equations. It must have 
determinant zero. This gives a polynomial equation in 
A, the eigenvalue. Write /(A) = det(A'( m ) (A)), which 
shows the explicit dependence on the eigenvalue A. We 
also have g(X) = det(A( m )(A)) = 0, the original charac- 
teristic equation. The coefficients in /(A) and g(X) are 
functions of the variables A^ . The necessary and suf- 
ficient condition that H has an eigenvector in a DFS is 
that /(A) and g(A) have a common root: that is, the re- 
sultant R(A) = 0. It can be shown that the resultant 
does not vanish identically. All possible Hamiltonians, 
parametrized by the real numbers A,j , constitute a man- 
ifold of dimension N(N + l)/2. The Hamiltonians which 
have an eigenvector in a DFS in V m lie on a submani- 
fold of dimension strictly less than N(N + l)/2. Hence, 
using Sard's theorem [1 51 ] and generalizing to Hermitian 
matrices we arrive at the following theorem. 

Theorem 1 Let X be the A^-dimensional real mani- 
fold of the parameters characterizing the possible Hamil- 
tonians in the Lindblad master equation. Let S be the 
set of values of the parameters which characterize Hamil- 
tonians that have at least one DF state other than the 
ground state evolving into a DFS state at all times. Then 
S is of measure zero in X. 

Proof Let S m be those members S that correspond to 
Hamiltonians that have some DF states in a fixed weight 
space V m , m > evolving into DF states at all times. 
From Proposition 1, we see that this is equivalent to the 
condition that the Hamiltonians in S m have an eigenvec- 
tor in the DFS in V m . From the preceding discussion, 
S m has measure zero in X. Since S = U m S m , S too has 
measure zero. 

Informally, we can say that for almost all Hamiltonians 
there is no DF state other than the ground state which 
evolves into other DF states at all times. Therefore, the 
best we can hope for is to seek states which remain DF 
up to some orders of perturbation. 
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Corollary to Proposition 1 implies that to get DF 
states we should look for states that are annihilated 
by operators S_, [H, S_], [H[H, S_]], .... We have seen 
that in general DFSs are not invariant under H. So, 
we seek invariance up to certain orders. The condi- 
tion for states \a) that are invariant up to first order 
(in H) is that S_ \a) = [H,S_]|a) = 0. Similarly, 
the second order condition is S_ \a) = [H, S_] \a) = 
[H, [H, S_]] | a) = 0. We assume that the initial 
state p(0) = |^} (ip\ is a DF state. Then p(t) = 
e -it{adH+iL D ) p(Qj anc j jf tne evolution is unitary, p'(t) = 
e~ ltIi \ip) e ItH . Hence, we take F 2 (p' (t) , p(t)) = 
(V>|e itH e _it ^ H+iir, V(0)e- itH \yj)a,s a measure of de- 
viation from unitary evolution [16(. Writing adH = 
X, we have F 2 = (M e ^x e -it(x+iL D ) Using tnc 

Zassenhaus formula [L7J , the (exact) fidelity can be writ- 
ten F 2 = (^|e L "e- i * 2 [ X ' L "]/ 2 e-* 3 ^(^ L «'---p(0)|V'). 
where Z 3 (X,L D ) = [X, [X, L D ]]/6 + i[L D , [L D , X]]/3) 
We consider the following three approximate cases. 

Case 1 (weak unitary part). The Hamiltonian 
H is replaced by eH, and treated as a pertur- 
bation. Then up to first order in e, F 2 = 

e tLo e -it 2 [L D ,X]/2 . . . e c fc t*(BdL D )»- 1 (JQ p ( ) |^ where 

Cfc = {-l) k - 1 e{k-l)i/k\, k > 1. The state p(Q) = \ip) (ip\ 
is DF, so F 2 = + £e fe i fc (adZ fl ) fc - 1 (*)M0)[# 

For DF-states (V>| (ad J L D ) fe - 1 (X)p(0) |?/>) = and i* 12 = 

i SGI. 

Case £ (strong unitary part). The Lindblad 

part is treated perturbation to give p = 

(— iadH + eLrj)p- Then up to first order in 

e, F 2 = ^\ e tL D ^[X.L D }/2... e d k t\ a dX)"-\L) p{ Q ) ^ ) 

where dfc = (~i) k e/k\, k > 1. Computing the se- 
ries expansions to first order, we get F 2 = (ip\ (1 + 

dkt k (a,&X) k ~ 1 {L D ))p{0) \ip). Unlike the previous case 
Ld(p(0)) = does not guarantee F 2 = 1 for all times. 
There are two possible courses of action. First, if the 
time scales of H are much shorter than Lp, zeroth or- 
der may suffice. Thus, instead of talking the initial state 
p(0) as an eigenstate of Lindblad generators, we use an 
eigenstate of H. This method is used in Rcf. [l||. 

Case 3 (short-time expansions). Note that the Zassen- 
haus formula shows that DFSs are stable up to first or- 
der in t. If we require stability up to order t k , k > 
1, we can work in the smaller subspaces W satisfying 
adH m (S,) • W = 0, m = 0, . . . , k - 1 V i. Since we are 
dealing with finite-dimensional Hubert spaces, it suffices 
to require stability up to order n, the dimension of the 
system space. However, we have seen that this is not pos- 
sible in general (Theorem 1), so we have to be satisfied 
with smaller k. 

We consider the universality of subspaces for external 
Hamiltonians that implement quantum gates. Let H c de- 
note the external Hamiltonian and H<; = Af 7 -S! ) _S'L. 
The equation of motion is p = (— i(adH c + adH^) + 
Ld)p = (Y + X + Lrj)p. Since we aim to work in DFSs, 



they must be invariant under any external Hamiltonian. 
Thus, [H c , Hi] = 0, and [Y,L D ] = 0. Assuming that 
the characteristic time-scale of H c is much smaller than 
other operators (essential for feasible computation) , then 
using the short time expansion we can construct robust 
subspaces for computation. For example, if we wish to 
terminate at third order, then a subspace W satisfying 
L D ■ W = [H d ,L d ] ■ W = [Hd, \H d ,L D ]] ■ W = will 
suffice. 

We have shown that in Markovian systems the Lamb- 
shift-type system Hamiltonian generally causes transi- 
tions from DF states to non-DF states. The results pre- 
sented here emphasize the importance of accounting for 
both unitary and nonunitary evolution in passive quan- 
tum error correction. Note that the results presented 
here, as long as the Lindblad operators do not cause 
transitions between irreps, can be extended to the finite 
temperature case. 
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